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Finite magnesium oxide nanotubes are investigated. Stacks of four parallel squares, hexagons, 
octagons, and decagons are constructed and studied by the pseudopotential density functional the- 
ory within the local-density approximation. Optimized structures are slightly distorted stacks of 
(N ' poly gons. These clusters are insulators and the band gap of 8.5 eV is constant over an investigated 

range of the diameters of stacked polygonal rings. Using the Lowdin population analysis a charge 
transfer towards the oxygen atoms is estimated as 1.4, which indicates that the mixed ionocovalent 
bonding exists in investigated MgO nanotubes. 

^ \ PACS numbers: 61.46.-fw, 73.22.-f; 36.40.Cg, 81.07.De 
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I. INTRODUCTION 



Studies of nanotubes and nanowires are important for recent research in nanotechnology. Experimental and theo- 
retical investigations of cylindrical nanostructures for many different organic and inorganic materials already existi 
As new arrangements of atoms these structures are interesting for basic science. There is also a possibility to use cylin- 
drical nanostructures as electrical and mechanical parts in nanodevices. Some technological applications of nanotubes 
and nanowires already exist, most notably for carbon nanotubes. Magnesium oxide is important in applications. 
It is known as an inert material with a high melting point, as a catalyst in important chemical reactions, and as a 
good substrate for many chemicals, for example metals, group III-V elements, and high-temperature superconductors. 
MgO is also used as an optical material and a component of composites. A bulk magnesium oxide crystallizes in the 
rocksalt structure. It is an insulating ionic simple oxide. 

Properties of materials at a finite scale of lengths are often different in comparison to their bulk. For example, 
a pronounced covalent contribution to the ionic bonding exists in small MgO nanoparticles, whereas almost pure 
l— — '■ ionic bonds are typical for the bulk of this compound. Therefore, a crossover from a mixed covalent and ionic 
I ' bonding to dominantly ionic one occurs in MgO nanostructures of an intermediate size. Small clusters of MgO have 
, been investigated by several experimental 2 ^ and theoretical me thodsi2i2i£i£iL2i2ii2iiiii2iI£ A study of MgO clusters 
f***. ' and their nucleation from the gas phase has been used to explore one possible process of a cosmic iron-magnesium- 
silicate dust formation in the last phase of life for some stars^ Mass spectroscopy and collision-induced-fragmentation 
measurements have been performed on sputtered MgO cluster ions. 2 Results have been shown that for small (MgO)„ 
clusters especially pronounced peaks exist at intervals of n = 3 (i.e., for n = 6,9, 12, and 15). These peaks indicate 
the most stable structures. Laser-ionization time-of-flight mass spectrometry measurements on (MgO)„ clusters 
produced in a gas aggregation source have also been carried out£ The clusters of enhanced stabilities were found 
for n — 2,4,6,9,12, and 15. Ziemann and Castleman have been performed calculations using two polarizable ion 
shell models^ Their results have been shown that the most stable structures are compact cubelike ones in the 
model with ion charges ±1. Conversely, spherical geometries composed of hexagons and squares have been found 
I ■ in the second model with ion charges ±2. Hollow and compact cubelike clusters have been also found in various 
other calculations4iiii§iL2i2ii2iiiii2iiii Cubelike structures are similar to pieces of the bulk MgO face-centered cubic 
lattice. These cubelike clusters prevail when the number of particles increases. Planar geometries, as well as spherical 
and cylindrical three-dimensional hollow ones, have also been found. A selfconsistent tight binding calculations of 
structural and electronic properties of (MgO)„ clusters have been carried out by Moukouri and Noguera*^ They have 
been calculated that planar polygonal rings are stable configurations up to n = 6. The (MgO)„ clusters of the same 
size have been investigated by the density functional DMOL package.^ These results have been shown that investigated 
5h ' small (MgO)„ clusters are the most stable for n = 3 and n = 6. Structures of MgO clusters have been studied by a 
genetic algorithm optimization method within a simple empirical rigid ion model^. The variation of cluster geometry 
with the charges ±q on the ions have been investigated. When the charge changes from q = 1 to q = 2, the structure 
evolves from compact cubelike to hollow pseudo-spherical one. Stacks of polygonal (MgO)„ rings have been found 
by the Hartree-Fock calculations for small clusterSf2ii£ and by optimization using the Born-Mayer potential^ Using 
an empirical compressible-ion interaction model potential and coordination number-dependent oxide polarizabilities 
Wilson has been found stable (MgO)„ nanotubes for n < 30iii Electronic and related structural properties of these 
MgO cylindrical hollow structures, i.e., nanotube clusters, have not been investigated. 

Several methods of fabricating MgO nanowires have been developed. Boron oxide assisted catalytic growth has 
been used to synthesize MgO nanowires with uniform diameter distribution of about 20 nm, and a typical length of 
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several tens of micrometers ii* Sometimes, nanowires with length of several hundred micrometers have been observed. 
X-ray diffraction, scanning electron microscopy, and high-resolution transmission electron microscopy techniques have 
been revealed that these nanowires are single-crystal fee structures. MgO nanowires with diameters of (15 — 20) nm, 
lengths of up to several hundred micrometers, and the fee structure have been also prepared using a vapor-phase 
precursor method starting from the MgB2 powders^ Cubic MgO nanowires with diameters of about 100 nm and 
lengths of up to several micrometers have been synthesized by heating the Mg powder above its melting point in a 
flow of 20% O2 and 80% Arii& Recently MgO nanotubes have been synthesized using a liquid metal assisted growthiii 
They have been filled with gallium and this nanoscopic device is proposed as a thermometer with better properties 
than previously constructed one based on a carbon nanotube>i& The fabricated MgO nanotubes are cubic, several 
micrometers long and exhibit both interior and exterior square cross sections. Their outer diameters are (30 — 100) 
nm, whereas interior ones are (20 — 60) nm. MgO nanorods have been used to produce columnar defects in high- 
temperature superconductors^^ They could pin magnetic flux lines and prevent thermally activated flux flow. Small 
MgO nanotubes studied in this work are also candidates for applications as components in the bulk, film, or wire 
composites of superconducting and other materials. 

In this work MgO nanotubes having between 16 and 40 atoms, and made of stacks of four polygons are investigated. 
By stacking four polygons on top of each other and optimizing these structures, slightly distorted nanotubes of the 
length from 0.5782 nm to 0.5925 nm are formed. The structure, energetics, and electronic properties of these MgO 
nanotube clusters are studied using computational methods based on a pseudopotential density functional theory in 
the local density approximation (LDA). The rest of the paper is organized as follows. Section ITT1 gives an outline of 
the method. In Sec. Illll thc results and discussion are presented. Conclusions are summarized in Sec. II VI 

II. COMPUTATIONAL METHOD 

In this work ab initio density functional theory calculations^*^ are performed to study the structural and elec- 
tronic properties of finite magnesium oxide nanotubes. The plane-wave pseudopotential method is applied^iS4 The 
density functional theory within the same plane-wave pseudopotential framework has already been used to calculate 
properties of the bulk MgO under high-temperature and high-pressure conditions which are important for geophysics 
and planetary physicsi2£*2& Recent density functional theory study of the high-pressure behavior of MgO at the gener- 
alized gradient approximation level (GGA)22i has shown an agreement with the LDA calculations of Refsi2£i2& Within 
the same plane-wave pseudopotential model as used in this work, it has been shown that the LDA approximation 
for the bulk of MgO and other alkaline-earth oxides yields good agreement with available structural and vibrational 
experimental data^ 

In this calculation LDA and parametrization of Perdew and Zunger for the exchange-correlation energy is applied^ 
The pseudopotcntials generated by the method of von Barth and Car for magnesiumSiiSi, and by one of Troullicr 
and Martins for oxygen^i are used. The configurations 2s 2 2p 4 of O, and 3s 2 of Mg are taken as valence states. 
The calculations are performed with a kinetic-energy cutoff of 60 Ry. In the studies of MgO under geophysical 
conditional*^ the cutoff of 90 Ry has been used, and this is important for a behavior of materials at high pressures. 
In the studies presented in this work (as well as in Refi 2 ^) it has been found that the results converge for the kinetic- 
energy cutoffs below 60 Ry. The cluster is positioned in the center of a large supercell of side 30 a.u. Such a choice 
ensures that the minimum distance between two clusters in the neighboring cells is in the worst case larger than 9.375 
A, which is found to be sufficient to prevent their interaction. The Brillouin zone is sampled using the T-point. 

Structural relaxation for MgO nanotube clusters is carried out by performing a series of self-consistent calculations 
and computing the forces on atoms. Rings and stable stacks of rings for MgO clusters have already been suggested by 
empirical and Hartree-Fock calculations £££2*^12^ Therefore, here a cylindrical distribution of atoms in the form of 
stacks of polygons is assumed and interatomic distances are varied to find a minimum energy configuration. Starting 
nanotubes are constructed as strictly parallel stacks of regular polygons of the sides of 0.2 nm, and at the distance 
of 0.2 nm. The experimental nearest-neighbor distance in the bulk magnesium oxide is 0.21056 nm. 32 Calculations 
have been shown that this distance decreases in nanoparticlesi 3 ^^*SiLSi2ii2iii4i2ii^ The experimental distance between 
atoms in the dimer MgO is 0.1749 nmS Stacks of four parallel squares (4 x 4), hexagons (6 x 4), octagons (8 x 4), 
and decagons (10 x 4) are studied. Therefore, the initial lengths of nanotubes are 0.6 nm. Several similar nanotubes 
are made for each size by breaking the symmetry of an ideal starting configuration by small random displacements 
of atoms. This procedure yields negligible difference in the final energy minima. The geometry optimizations are 
performed using the Broyden-Fletcher-Goldfarb-Shanno method of minimization. All atoms are allowed to relax 
without any imposed constraint. The optimizations are carried out till the forces on the atoms are below 10 -4 a.u. 
For geometries far from the minimum, the self-consistency cycle continues till the energy differences are less than 10~ 8 
a.u. Close to the minimum smaller self-consistency thresholds are taken (down to 10 -10 a.u.) in order to produce 
correct forces and small energy differences involved in the relaxation^ MgO nanowires and nanotubes synthesized in 
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TABLE I: The geometry of optimized nanotubes, a is the length of the polygon side, I is the distance between neighboring 
polygonal rings, L is the length of the tube (i.e., the distance between opposite atoms on the top and bottom polygons). Indices 
o and i respectively label outer (i.e., when at least one atom is in the top and bottom polygons) and interior bonds. Units are 
nm. The range of values is given if more than two values appear for a certain distance. 



Nanotube 


4x4 


6x4 


8x4 


10 x 4 


a 


0.1937 


0.1895; 0.1900 


0.1883; 0.1884 


0.1874-0.1880 


di 


0.2101 


0.2029 


0.2007; 0.2009 


0.1997-0.2004 


la 


0.1919; 0.1956 


0.1950-0.2000 


0.1961; 0.2017 


0.1963-0.2019 


u 


0.1922 


0.1954 


0.1961 


0.1969 


L 


0.5782 


0.5883; 0.5886 


0.5918 


0.5925 


TABLE II: 


Total energies per MgO unit. 


D is the diameter of the enclosed 


cylinder for an initial regular structure. 


Nanotube 




Diameter (nm) 




Total energies (eV) 


4x4 




0.282845 




-464.6981 


6x4 




0.400000 




-465.1345 


8x4 




0.522603 




-465.2519 


10 x 4 




0.647249 




-465.2947 



T^ e fg | i4 i i5 4 i6 i i7 are typi ca iiy several tens of micrometers long, and their outer diameter is (15 — 100) nm. These sizes 
are not accessible for computers used presently. It is found that infinite periodic MgO tubes in the form of stacks of 
polygons similar to nanotube clusters studied here, but in the superlattice geometry, are not stable. 

III. RESULTS AND DISCUSSION 

Figure Q presents the structures of MgO nanotubes which correspond to the minimum energy. This visualization 
is performed using the Rasmol packagei^S The results show that stacked ring structures are stable during a pseu- 
dopotential density functional theory minimization. As presented in Fig. ^ the bonds and angles slightly change 
during the geometry optimization. The deviation from planarity of initial polygons is small. Similar, but slightly 
more distorted structures have been obtained by Wilson using an empirical compressible ion potential model^i 

It is found here that the square cross-section of the smallest 4x4 nanotube evolves to the rhombus. For the top 
and bottom rhombuses the O-Mg-0 angle is 96.3°, whereas the Mg-O-Mg angle is 83.5°. The 6 x 4, 8 x 4, and 10 x 4 
nanotubes after optimizations have a barrel shape as shown in Fig. ^ The optimized structural parameters are shown 
in Table. [I] The length of nanotubes increases with the number of atoms in the cluster. In optimized structures 
majority of the Mg-0 distances decrease in comparison with initial ones of 0.2 nm. Total energies as a function of the 
initial diameter of MgO cylinders are shown in Table. ^ 

The electronic density of states (DOS) for the 8x4 nanotube is shown in Fig. [21 The separation of bands is about 
8.5 eV for all investigated nanotubes. When the diameter of a nanotube increases DOS does not change substantially, 
and the peaks only increase in their heights. A gap of 5.5 eV independent of the tube diameter and chirality has been, 
for example, found for tubular boron nitride nanotubes^ The band gap value of about 9 eV has been calculated 
for the 64 atoms cubelike MgO cluster using the density functional DMol packaged However, there the decrease 
of the band separation when the size of cubes increases has been reported. It has been calculated that for the 
(MgO)4 cubelike cluster the band separation is 12 eV. The increase of the number of atoms changes the number of 
neighbors for a cube cluster geometry where the coordination in larger clusters approach the bulk value. Conversely, 
the coordination in the (MgO)„ nanotube clusters studied here does not approach the bulk value for large n. A strong 
shape dependence of electronic properties in small MgO clusters has been found by Moukouri and Noguera using a 
tight binding calculationi^ It is known that both LDA and GGA exchange correlation functionals underestimate the 
electronic band gap. The band gap of 7.7 eV has been calculated in the GW approximation to the bulk electronic 
structure of MgO, in agreement with the experimental value of 7.8 eVml The LDA approximation used within various 
methods of electronic structure calculations gives (4.19 — 6.6) eV for the bulk band gap of MgO. 37 Therefore, the 
band gaps of 8.5 eV calculated in this work for MgO nanotubes within LDA, and 9 eV found in Refii£ for the cubelike 
(MgO)32 cluster in the same approximation, are much larger than the values for the LDA bulk gap reported in the 
literature. An increase of the band gap in nanoparticles in comparison to the bulk occurs also for other materials. 
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FIG. 1: The optimized geometries of MgO nanotubes: (a) 4 x 4, (b) 6 x 4, (c) 8x4, (d) 10 x 4. Dark and light balls are used 
to represent the Mg and O ions, respectively. Side views and views along the tube axis are shown. 
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FIG. 2: Density of states for the 8x4 MgO nanotube. DOS for other nanotubes is similar. 
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TABLE III: The results of the Lowdin population analysis. Average charge transfer SQ to oxygen atoms is shown for all 
nanotubes. Q is the charge on the atom whose number of neighbors is Z (initial valence charge is 2 for Mg, and 6 for O). S is 
the spilling parameter. 



Nanotube 4x4 6x4 8x4 10 x 4 

SQ 1.39 1.40 1.41 1.42 

Q(Mg, Z=3) 0.44 0.42 0.40 0.39 

Q(Mg, Z=4) 0.44 0.43 0.42 0.41 

Q(0, Z=3) 7.37 7.39 7.40 7.41 

Q(0, Z=4) 7.40 7.41 7.42 7.43 

S 0.02 0.02 0.02 0.02 



The origin of this phenomenon is a confinement in finite systems which increases the separation between energy levels. 
The value of the band gap of various MgO clusters has not been measured, and this subject certainly deserves an 
investigation. 

The Lowdin population analysis has been carried out^& The charges on the atoms and the average charge transfers 
are shown in Table. ILTT1 The spilling parameter S measures the difference between the plane- wave eigenstates and their 
projection into the atomic basis, i.e., the amount of the total charge lost in the projecting process^ Identical atoms 
with different coordination on the same nanotube sometimes have a different charge. The large charge transfer, but 
still less than two, shows that the mixed ionic and covalent bonding exist in these nanotube clusters. As in RefsASil, 
the term ionocovalent is used here to label the mixed ionic and covalent bonding. Although in a population analysis 
and other calculations involving the charge of ions the values depend upon the calculation model and method, a short 
discussion of several results follows. This point strengthens a connection between structural features of nanoparticles 
and the charges on the ions. The analysis of cubelike MgO clusters^ has been shown that charges depend on the size 
of nanoparticles and on the coordination number of atoms. Moukouri and Noguera have been found a large charge 
reduction as the coordination number decreases yielding to the charge in the MgO dimer of only half that of the 
bulk ions££ In an early study based on the rigid ion approximation and two polarizablc ion shell models, Zicmann 
and Castleman have found that cubelike structures of (MgO) n correspond to the model with ion charges ±li£ In 
their second model with ion charges ±2 spherical geometries composed of hexagons and squares have been optimized. 
Similar results have been obtained by a genetic algorithm method. 8 Using an empirical potential with fluctuating 
charges Calvo has been found that for cubelike (MgO) n clusters, when n changes from small values to the bulk 
regime, the charge transfer increases from 1 to 2^ The crossover from ionocovalent to ionic bonds has been estimated 
to occur in these clusters at (300 ± 100) MgO molecules. In the nanotube geometry studied here the coordination 
numbers are the same for all clusters and for larger clusters do not approach the bulk value. Therefore, the charge 
transfer of about 1.4, and the ionocovalent bonds are typical for these MgO nanotube clusters. The contour graphical 
representation of the charge density for the 6x4 cluster is presented in Fig. [3] This visualization is performed using 
the XCrySDen package. 40 Charges are mostly located on the oxygen atoms, but the covalent contribution to the ionic 
bonding exists. 

IV. CONCLUSIONS 

MgO nanotubes consisting of four stacked planar n-polygonal rings (n — 4, 6, 8, 10) are investigated by the first- 
principles density functional theory method in the pseudopotential approximation. These calculations show that MgO 
nanotube clusters are stable slightly distorted stacks of polygons. They have a large insulating band gap of 8.5 eV 
independently of their radius. The Lowdin population analysis indicates that the average charge transfer to oxygen 
atoms is 1.4. Therefore, the covalent contribution to the ionic bonding exist in these MgO nanotubes. 

Experimental studies of various MgO clusters, for example done by photoemission and electron microscopy tech- 
niques, are desirable. Recently an atomic resolution imaging of the Mg(001) surface has been reached using a dynamic 
scanning force microscopy^i This experiment shows that similar imaging of MgO nanotubes is also feasible. It will 
be interesting to prepare, study, and use multiwalled MgO nanowires similar to carboni and gold 4 ^ ones. MgO 
nanotubes are cylindrical structures with large interstitial regions made of a chemically inert material. Therefore, 
they are promising for use as nanocapsules. Submicrometer-sized MgO cubes have already been used as a template 
for growth of carbon nanotubes and these structures were imaged by electron microscopy^ It is possible to prepare 
similar templates based on MgO nanotubes. Rectangular MgO nanorods have been recently used for the growth of 
embedded In nanowires^ Such MgO composite nanostructures are promising candidates for antioxidation coating of 



FIG. 3: The charge density contour plot on a plane perpendicular to the axis of the 6x4 cluster and passing through the upper 
hexagon. The 0.01 eja% isovalue is shown. The small dark spheres are O, whereas large light ones are Mg atoms. 



metal nanowires, and as insulating layers in nanoelectronic devices. As for the BN nanotubes^ 6 ., a radius independent 
band gap calculated in this work offers a possibility to use grown samples of MgO nanotubes containing different sizes, 
but with similar electronic properties. Therefore, many possibilities for synthesis and experimental investigations, as 
well as for applications in various fields of nanotechnology exist. Some important applications of MgO cylindrical 
nanostructures have already been describedii2ii2iS£ 
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